Vortex solid phase with frozen undulations in superconducting Josephson-junction 

arrays in external magnetic fields 
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A vortex solid with self-generated randomness is found theoretically in a frustrated Josephson 
junction array (JJA) under external magnetic field with anisotropic couplings. Vorticies induced by 
external magnetic field develop stripes parallel to the direction of weaker coupling. It is shown ana- 
lytically that there is a continuous, gapless band of metastable states in which stripes are deformed 
randomly by transverse undulation. The vortex solid with the frozen undulation in a metastable 
state freely slides along the direction of stronger coupling, thereby destroying ordering of phases 
even at zero temperature, but is jammed along the direction of weaker coupling. 

PACS numbers: 61.43.Fs,62.20.Qp,74.81.Fa,74.25.Qt 



Solids are systems with rigidity ranging from crystals, 
uasi-crystals to glasses and jammed granular matters 
l| . Generally it becomes more challenging to understand 
mechanism of the formations of less periodic solids distin- 
guishing them from liquids Q • A useful guiding concept 
to study non-crystalline solids is frustration which in- 
hibits simple crystallization. Imagine that there is a trick 
to "inject" dislocations artificially into a crystal from out- 
side and that their density can be controlled at will. Such 
a system will provide a very interesting ground to study 
consequences of geometrical frustration, especially real- 
ization of self- generated randomness or glassiness without 
quenched disorder 0,3- Quite remarkably the Joseph- 
son junction array (JJA) under external magnetic field 
realizes such an ideal situation . Furthermore, trans- 
port properties of the JJA under external current can be 
regarded as "rheology" under external shear ljj 

JJA is a regular network of superconducting islands 
connected with each other by Josephson junction in the 
form, say, of a square lattice of size N = L x L as shown 
m Fig. CD a) [lEf. The phases 9i of the superconducting 
order parameter on the islands i = 1,2, ... ,N interact 
with each other via Josephson coupling. Under magnetic 
field B, which can be varied at will in experiments, the 
number density / = Ba 2 /(j>Q of vorticies (dislocations) 
can be forced into the configuration of the phases. Here 



a is the area of a plaquette and <po is the flux quantum. 

An interesting connection of JJA to the problem of 
frictions provides valuable insights. The frustrated JJA 
becomes essentially equivalent to the Frenkel-Kontorova 
(FK) model [11] and the two-chain model of Matsukawa 
and Fukuyama (MF) [l3| in one dimensional limit (TV 



L X 2), i. e. on the ladder lattice [ll| 14]. These one- 
dimensional systems are known to undergo a kind of jam- 
ming or frictional transition at zero temperature T = 0, 
known as the Aubry's transition 15, 16[ at a critical value 



of the strength of coupling A between two surfaces which 
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FIG. 1: Vortex patterns in an irrationally frustrated JJA un- 
der external magnetic field with anisotropic coupling. Here 
A = 1.5 so that the coupling is stronger along y direction. 
Such an anisotropic JJA can be fabricated in laboratory by 
lithography technique(s) [jj. Panel a) displays the JJA on a 
square lattice. A fraction / = 21/55, which approximates an 
irrational number (3 — v5)/2 = 0.381966.., of the plaquettes 
are occupied by vorticies with charge 1 — / represented by 
filled squares. Panel b) displays an equilibrium vortex pat- 
tern at T = 0.2 and Panel c) displays that at a nearby energy 
minimum reached via an energy descent algorithm. 



are incommensurate with respect to each other. Then one 
would naturally be led to consider irrationally frustrated 
anisotropic JJA (lp| with 1) irrational vortex density / 
and 2) anisotropic couplings into x and y directions, 
say 1 into x direction and A into y direction. 

In this Letter we study the ground state as well as 
low-lying states of the irrationally frustrated JJA with 
sufficiently strong anisotropy A ^> 1. Numerically we 
found vortex stripes parallel to the direction of weaker 
coupling. In addition to the ground state in which vortex 
stripes are straight, we found numerous metastable states 
with different realizations of transverse undulation of the 
stripes as shown in Fig. [T] By a perturbative analysis in 
series of 1/A starting from infinite anisotropy limit A = 
oo, we are able to reproduce a family of such low- lying 
metastable states including the ground state analytically. 
The coexistence of sliding and jamming in the system 
fiol ] is proved from the analytically constructed ground 
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state and the low-lying states. Because of the sliding, the 
phases remain disordered even at T = for irrational /, 
in sharp contrast to JJA with rational / where not only 
vorticies (chiralities for / = 1/2) but also phases exhibit 
(quasi-)long ranged order at T > (l7j . 

Model To simplify notations we label the vertices (su- 
perconducting islands) as i = 1, 2, . . . , N whose position 
in the real space is given by (n.j,mj). The properties of 
the JJA under the transverse magnetic field are known 
to be described by an effective classical Hamiltonian [f|, 



H 



cos(^)- 

<t,i>||x-axis 



A 



cos(^) (1) 



<iJ>||y-axlS 



where < i,j > denotes nearest neighbor and tpij the 
gauge-invariant phase difference, ipij = 6i — 0j — Aij. 
The temperature T is defined in a unit with fee = 1. For 
the anisotropy A, we need to consider only A > 1 by sym- 
metry. The vector potential Ay (= —Aji) is defined such 
that directed sum of them around each plaquette is 2nf. 

Vortex charge v t of the vortex at the plaquette asso- 
ciated with the i-th vertex is defined by taking directed 
sum of (ijjij /2tt — [ipij/2ir} n ) on the junctions around the 
plaquette. Here [x] n denotes the nearest integer of the 
real variable x. It takes values . . . , — 1 — /, — /, 1 — /, . . .. 
We use periodic boundary conditions so that the total 
vortex charge is enforced to be zero (charge neutrality). 

It has been proposed that superconducting glass may 
be realized if / is irrational [|| . While J JAs with rational 
/ develop periodic vortex lattices [T3, [HI ; such simple or- 
derings may be avoided with irrational /. Indeed equilib- 
rium relaxations were similar to the primary relaxations 
observed in typical fragile supercooled liquids [l9 . Such 
a system is called as irrationally frustrated JJA. 25f 

Numerical methods In numerical simulations, we 
used a series of rational numbers p/q = 5/13, 8/21, 
13/34, 21/55, 34/89, 55/144, 89/233 for the filling /, 
which approximate an irrational number / = (3 — 
v / 5)/2 = 0.38196601... Square lattices of size L x L 
with periodic boundary conditions for both directions are 
used. We choose L — q so that the ratio f — p/q con- 
verges to the target irrational number in L —¥ oo limit. 

To generate the equilibrium ensemble, we used ex- 
change Monte Carlo (MC) simulations combined with the 
over-relaxation method [20( performed on systems with 
L = 13 — 89 using 20 — 120 temperatures in the temper- 
ature range T = 0.2 - 0.4. We used 10 5 - 10 6 MC steps 
for the equilibration and observations. 

Stacked undulating vortex stripes As shown in 
FiglU the vorticies develop undulated stripes parallel to 
the direction of weaker coupling at low temperatures. 
The formation of the stripes is reasonable because the 
repulsive interactions between vorticies are anisotropic. 
A remarkable feature is that the stripes are stacked quite 
regularly along the stronger coupling as shown in Fig[TJ 
c) in the nearby energy minima obtained via a simple en- 
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FIG. 2: Structure factor of vorticies. a) displays the cross- 
sections of the structure factor S(q x ,q y ) (L = 55) with ther- 
mal average (thick lines) and at the energy minimum shown 
in Fig.[T]c) (thin lines), b) displays the amplitude of the peak 
of the structure factor S* = S(n,2nf) with thermal average 
(S t ), average over the energy minima [S*] m i n . and variance of 
the minima-to-minima fluctuation ^/[5»] m in. — [5*] min • Here 
the average over minima [. . .] m i n is taken over 100 energy- 
minima obtained by independent initial conditions. 



crgy descent algorithm. Starting from different thermal- 
ized configurations we obtained numerous energy min- 
ima similar to the one shown in FigfT] c) but with differ- 
ent realizations of the transverse undulation. The nearly 
perfect stacking of the stripes strongly suggests that the 
energy barrier in going from one to another realization 
of significantly different undulation of vorticies (disloca- 
tions), which necessarily involve large number of plastic 
deformations, diverges with the system size so that the 
ergordicity is broken. This feature is markedly differ- 
ent from usual undulations found, for example, in liquid 
crystals which are fluctuating dynamically [2l|. 

The stacked undulation is manifested in the struc- 
ture factor of the vorticies. As shown in Fig. [2J the 
structure factor S(q x ,q y ) exhibits prominent peaks at 
(lx,qy) = (tT) 27t/) and (tt, 2tt(1 — /)) whose height scales 
with the system size as N = L 2 as usual Bragg peaks 
do. However the profile of the peak is peculiar: it de- 
cays sharply along q y reflecting the stacking but decays 
slowly by a power law \ir — q x \~ 2 along q x reflecting the 
transverse undulation. Here we emphasize again that the 
transverse undulation is frozen in time. The frozen-in 
randomness is manifested in the minimum-to-minimum 
fluctuations of the structure factor shown in Figf2]a) and 
b). Note that the variance of the fluctuation as well as 
the average grows linearly with the system size N mean- 
ing that the structure factor is not self-averaging. 

Analytic construction of the ground state Let us 
now turn to explicit construction of low-lying states by 
an analytical approach. To this end we propose a non- 
trivial ansatz for the ground state using the notion of 
the so called hull functions developed in the studies of 
the FK and MF models fT^4l6l ] . We propose that the 
gauge invariant phase differences ipij across the Joseph- 
son junctions, in the low-lying states of the anisotropic 
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JJA (A > 1) can be represented as, 

ip(x,y)(x+i, v ) = 0x[y+a{x)} ip( x ,y)(x.y+i) = <f> y \y+a(x)] 

(2) 

Here (f> x [y] and 4> y [y] are functions defined on the "folded 
coordinate" [y] — fy — 'mt(fy) where int(x) is the floor 
function. The folded coordinate takes values limited in 
the range < \y] < 1. Such a function is called as a hull 
function [12t 1 1 5I | . Note that if / is irrational, which we 
always assume in this work, the vertices of the JJA uni- 
formly fill the entire range of the folded coordinate [y] in 
the limit N — > oo. Thus we can treat [y] as a continuous 
variable. Moreover one can then easily extract distribu- 
tion of the values of the phase differences from the hull 
function because of the uniform distribution of [y] over 
the support < [y] < 1. 

An obvious constraint is that the directed sum over 
ipij around each plaqucttc must be — 2irf. In addition, 
the Josephson currents must be conserved at each vertex 
(force balance condition) in each energy minimum. Thus 
the following two conditions should hold, 

<t>M + Mv + s ( x )] - Mv + 1] - Mv) = -2*7 (3) 

- sin cj) x [y] + sin <f> y [y] = - sin <f> x [y - 8(x)] + sin <j> y [y - 1] (4) 

where 8{x) = a(x) — a(x — 1). 

Now our task is to look for the hull functions <fi x [y], 
4> y [y] and phase differences S(x) which satisfy the con- 
ditions on the plaquettes Eq. ([3]) and vertexes Eq. (0}. 
We solve this problem by performing a 1/ A expansion 
lrl 2^] around the infinite anisotropy limit A = oo. 

In the infinite anisotropy limit A —> oo the weaker cou- 
plings can be neglected so that we easily find <j> x [y] = 
(2[y] - 1)tt + 0(1/A) and <j> y [y] = 0(1/ X) which trivially 
satisfies Eq. and Eq. (UJ) (with A = oo). As such, the 
phase difference S(x) is not fixed at this stage. 

Let us consider first the ground state assuming that 
phase difference is uniform, i. e. d(x) = S. Us- 
ing the above results in Eq. (01 we find 1/A correc- 
tion term of 4> y [y], which in turn allows us to find 
1/A correction term of 4> x [y] through Eq. (J3J. In this 
way we obtained analytic form of the hull functions up 
to 0(A" 3 ) as, <t> x [y) = (2fy - 1)tt + \^^s x (l,y) + 
l^(\a 3 (d)\ 2 Sx (3,y)^3\ ai (S)\ 2 Sx (l,y)) and ^[y] = 

^s y (l,y) + «^ (\a 3 (S)\s y (3,y) - 3\ ai (6)\ Sy (l,y)) 
with a n (S) = (1 - e- a ^ s )/{\ - e~ i2n 'f), s x (n,y) = 
sin(2mr/i/) and s y (n,y) = sm(2mrfy + Arg a n (S)). For 
irrational /, we find the energy is minimized by choosing 
5 = S* = 1/(2/) and obtain the ground state energy E g , 

= _i_Mnf_MT ( i _ \ ai{s *)\2) + o(\-s). 

The vortex configuration of the ground state is indeed the 
stripes like Fig.[T]but without the transverse undulation. 

Band of undulated metastable states Next 
let us construct the low-lying states with transverse 
undulation of stripes shown in Fig. [T] Somewhat 
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FIG. 3: Configuration of the gauge invariant phase differ- 
ences across Josephson junctions in energy minima. Here the 
anisotropy is A = 1.5. In panel a) original y axis is "folded" to 
[y] — fy — int(/y) . The symbols are data of the gauge invari- 
ant phase differences tp x — ip( x ,y)( x +i,y) (filled symbols) and 
i'v = i'(x,y)(x,y+i) (open symbols) across Josephson junctions 
parallel to x and y axis. Each data set consists of data points 
on a 'column' (x, 1), (x, 2), . . . , (x, L) at an arbitrary chosen 
x (1 < x < L) in an arbitrary chosen energy minima respec- 
tively. Each data set is shifted globally by some a so that 
different data sets collapse on top of each other. The lines are 
analytically obtained hull functions of the ground state. In 
panel b) the data on the same energy minimum (L = 55) are 
plotted against "folded" a;-axis. 

surprisingly, we can solve Eq. ([3]) and Eq. d3J 
with arbitrary 5(x), finding, 4> x [y\ = (2fy — l)ir + 
[ Cl [S(x)] cos(2^ fy) + C 2 [S(x)]s x (l, y))] + 

0(A" 2 ) and 4>y[y] = (s y (l,y - S(x) + 1/(2/)) + fl „(l,y))+ 

0(A~ 2 ) withCi[5(x)] = -sin(27r/d(a;+l))-sin(27rjFd(a;)) 
and C 2 [S(x)} =-2 + cos{2irf6 (x + 1)) + cos(2irf6(x)). 
As the result the energy becomes, for irrational f,-^- + 

£l£3ipJ! J2t=A 1 + C0S ( 27r f 6 ( x ))) + °( x ^ 4 )- lt is evident 
that there exists a gapless, continuous spectrum of low- 
lying states each of which is parametrized by a function 
S(x). Assuming \S(x) — l/(2/)| <C 1, we obtain an one 
dimensional elastic Hamiltonian with an unusual elastic 
constant which grows linearly with the system size L. 

Let us emphasize that the undulated states with the 
arbitrary 6(x) are ensured to satisfy the force balance 
condition Eq. Thus the system trapped in an un- 
dulated state cannot relax spontaneously down to the 
ground state since they are metastable: the undulation is 
distinct from "phonons" by which vorticies cannot move. 

Now the unusual structure factor of the vorticies 
S(q x , q y ) in Fig.fJlcan be understood as follows. The con- 
figuration of the ground state is a function of the folded 
coordinate [y] so that S(q x ,q y ) must have peaks along 
the q y -axis at q y = 2-nf and 27r(l — /). On the other 
hand the phase shift of the hull function 5* = 1/(2/) 
along x-direction meaning that S(q x , q y \ must have a sin- 
gle peak along the q x axis at q x = n, reflecting the hor- 
izontal stripes. The power law tail \q x — ir\~ 2 naturally 
follows from the effective one-dimensional clastic Hamil- 
tonian for the transverse undulation obtained above. 
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In Fig. [3] we plot the phase differences across the junc- 
tions in energy minima obtained numerically and com- 
pare them with the hull functions obtained analytically. 
For simplicity we show here the hull function for the 
ground state disregarding small differences due to the un- 
dulation. Our perturbative result grasps well the overall 
features. The agreement will be further improved by tak- 
ing into account higher order terms in the 1/A expansion. 

A remarkable consequence of the analytic hull function 
is that the undulated vortex solid can freely slide into the 
y direction: Given an energy minimum described by the 
hull functions <f> x [y] and (f> y [y] , a family of solutions with 
exactly the same energy can be obtained through the op- 
eration [y] — > [y + a] with varying phase shift parameter 
a. Consequently the phases must remain disordered even 
at T = 0. As shown in Fig. [3] a) the phase differences 
across junctions parallel to the x-axis take all possible 
values in the range — ir < 4> x [y\ < tt meaning that the 
system can be sheared indefinitely along the a;-axis which 
amount to unidirectional sliding of the undulated vortex 
solid into the y-direction with fixed pattern. In contrast, 
the plot against "folded" £-axis shown in Fig. [3]b) ex- 
hibits no hint of a single valued, regular hull function. 
Also note that the distribution of the 4> y does not span 
the entire range — 7r < <t> y < vr which is needed to allow 
shear along y-axis mcasning jamming along x-axis. 

It is instructive to compare the above results with the 
FK model. In the FK model the hull function is proved 
to be an analytic function in the sliding phase A < A c 
but becomes non-analytic in the jammed phase A > A c 
[HI]. In the anisotropic J J A, the sliding and jamming 
are dual in the sense that they are simultaneously taking 
place along different axes. Indeed in [l(| it was found nu- 
merically that the shear-modulus is zero/finite along the 
direction of weaker/stronger coupling at T = 0. Further- 
more it was suggested that the symmetric point A = 1 is 
a critical point at zero temperature A C (T = 0) — 1 similar 
to the jamming point in granular matters Indeed re- 
cent studies at finite temperatures suggest T C (X = 1) = 
H, 24 1 . On the other hand the growth of the peak height 
of the vortex structure factor with the system size L 
(Fig. [2]) suggests T C (A) > at least at A > 1. 

The analogy to the FK model suggests that the ana- 
lytic hull functions 4>x(y) [y] becomes non-analytic at the 
critical point A = A c = 1 and remains non-analytic for 
A < 1. (The other way around for 4>x(y)[ x X)- However we 
speculate that the 1/A expansion, which yields only ana- 
lytic hull functions, remains stable up to A = 1. Then cer- 
tain "non-perturbative" solution(s) of Eq. ^ and Eq. 
must emerge at weaker anisotropy and make level cross- 
ing^) with the horizontal stripe state. 

To conclude we found undulated vortex stripes in 
irrationally frustrated Josephson junction array with 
anisotropic Josephson coupling theoretically. It will be 
very interesting to study critical properties of the system 
closer to the symmetric point where the present pertur- 



bative approach should break down. 
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